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Abstract. We derive a field-theoretical representation for the moments of the 
eigenstates in the generalized Anderson model. The representation is exact and can 
be used for the Anderson model with generic non-random hopping elements in any 
dimensions. We apply this method to the simplex model, for which the hopping 
amplitude between any two lattice sites is the same, and find that the eigenstates 
43 . are localized at any strength of disorder. Our analytical predictions are in excellent 

agreement with the results of numerical simulations. 
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1. Introduction 

The phenomenon of Anderson localization remains one of the most important and 
actively studied problem in the field of disordered quantum systems. Despite a great 
progress in understanding of this phenomenon achieved in the last fifty years (see p] 
and references therein), the number of available models affordable for non-perturbative 
analytical treatment is very limited. Among them are the original one-dimensional 
Anderson model [2] and its continuous version [31 HJ E] , quasi-one-dimensional models 
and their random matrix counterparts IE] , as well as models on a tree-like graphs 
Ell ttB 1121 [131 E|- A common feature of all such models is that they can be 
treated recursively. In the discrete case that means that it is possible to establish a 
relatively simple relation between properties of a system of size N and of size N + 1. 
A necessary condition for existence of such recursion relations is the absences of loops 
in the configuration space. Indeed, the existence of just a single loop, like in the case 
of the one-dimensional Anderson model with periodic boundary conditions, makes the 
recursive method inapplicable. 

In this paper we study the Anderson localization in a system containing a maximum 
number of loops for a given number of sites. This is the Anderson model on a <i-simplex, 
which is a generalization of the notion of a triangle in the two-dimensional space to the 
d- dimensional space. The Hamiltonian of the simplex model is similar to the standard 
Anderson model and it consits of the hopping term and the on-site random potential: 

H = T + V } T ij = ^ 1 V ij = v i 5 ij , i,j = l,...,N, (1) 

where N = d+ 1 is the total number of sites, 5ij is the Kronecker delta symbol and Vi are 
independent Gaussian distributed random variables with the zero mean value {v^ = 
and the variance (vf) = w 2 . The hopping amplitude between any two sites is equal to 
1/N and we set Tu = 1/N for the sake of convenience. 

In the absence of disorder the spectrum of the perfect system can be easily found 
if we notice that 

Tf = f , / = (1,1,...,1) T , (2) 

so that the vector fo, whose all componets are equal to one, is an eigenvector of T with 
the eigenvalue A = 1. On the other hand, any vector / orthogonal to /q satisfies the 
equation 

Tf = 0. (3) 

There exist N — 1 linearly independent vectors, which are orthogonal to fo, and hence 
we conclude that the second eigenvalue of T, A = 0, is (N — l)-fold degenerate. 

Thus the simplex model belongs to the class of the Anderson tight-binding models, 
whose clean analogs have flat bands of highly degenerate states. It is known that in 
the presence of disorder the eigenstates of such systems may become critical at weak 
disorder [15] or exhibit the Anderson transition at strong disorder [16] . For the simplex 
model one would expect that weak disorder should lift the degeneracy and facilitate the 
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emergence of extended eigenstates, as the hopping matrix connects any two sites of the 
lattice. However, we show that, in contrast to this expectation, all eigenstates, which 
were initially degenerate, are localized for any disorder strength in the thermodynamic 
limit N — > oo. To this end we calculate analytically the moments of the eigenstates as 
an explicit function of w and demonstrate that they remain finite in the limit iV — > oo. 

Another aim of this work is to present a new approach for calculation of the 
moments of the eigenstates suitable for Anderson models with non-trivial connectivity. 
Representing the disorder averaged moments of the eigenstates by the supersymmetric 
functional integral at the first step, we show how the most of the integration variables 
can be integrated out leading to the functional integral containing only a single real 
variable associated with every lattice site. As the constructed representation is exact 
and valid for the Anderson model with a generic hopping matrix, it can serve as a 
starting point for further investigations. In particular, for the simplex model it can be 
reduced just to a two-fold integral, which can be analyzed in the limit iV — > oo by the 
standard methods. 

The paper is organized as follows. In Section [2J we derive a general functional 
integral representation for the moments of the eigenstates of the Anderson model with 
a generic hopping matrix. In Section [3l we apply this representation to the simplex 
model and show how in this case it can be reduced to a two-fold integral. The latter 
is then evaluated in the limit N — > oo and compared with the results of numerical 
simulations in Section HI Finally we conclude the paper with a discussion of the origin 
of the localization in the simplex model and its generalizations in Section [5j 



2. Moments of the eigenstates 

The local moments I q {n) of the eigenstates on a given lattice site n are defined as follows 

I q (n) = TT^y E O^WI 2 " S(E -E a )), (4) 

where f a is an eigenstate of H corresponding to an eigenvalue E a , p(E) is the averaged 
density of states and (...) stands for the disorder averaging. 

The moments of the eigenstates defined above can be extracted from the singular 
part of a product of the diagonal elements of the retarded G R (E) = (E + it — H)^ 1 and 
the advanced G A (E) = (E — it — H)^ 1 Green's functions. Indeed, if we define K Lm (n,t) 
as 

^, m (n,6) = (Gi)'(Gl) m , /,m = l,2,..., (5) 
then the moments I q (n) can be found by taking the limit e — > [T7] 

'•<"> = C(W+"-2)'' SS <*<-<»■ £ » • « = ' + m - «o 

It is well known that the products of the Green's functions can be conveniently 
represented by the supersymmetric functional integrals [18] . In the case of a discrete 
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lattice of size N we introduce N supervectors associated with each lattice site: 

/ s R {i) \ 

XR{i) 

S A {i) 

\ Xa(i) ) 

whose components are two complex s R , sa and two Grassmann \Ri Xa variables. The 
product of the diagonal elements of the Green's functions can be written now as the 
Gaussian integral over supervectors $j [13, HE] : 

JV 



Ki tJn (n,e) 



l\ ml 



cxp 



p,q=l 

where $ f = (s* R , x*r, s* a , d$d& = -^2d' 2 s R d 2 s A dx R dx R dxAdx% the diagonal 
matrix A = diag(l, 1, — 1, — 1) and we define the scalar product of two supervectors 
as 

{%, = s* R {p)s R (q) + X R (p)Xn(q) ~ s A (p)s A (q) + X A {p)XA{q). (9) 

This representation enables us to perform averaging over the random diagonal part of 
the Hamiltonian ([I]) explicitly. Then the action of the functional integral (jSJ) takes the 
form 



J] d%d^ p (s* R (n)s R (n)) l (s* A (n)s A (n)r 
P =i 

N N 

i (H m - E5 pq )(%, %) - e YX* P , 

p=i 

1 ^2„_J2, 



N / 2 \ N 

p=l ^ ' p,q=l 



(10) 



It is convenient to introduce the generating function y($ n ,$JJ, which is obtained by 
integrating exp(— S[Q>, over all $ p except $ n : 



The lower and the upper limits for p is not written explicitly in order to lighten the 
notation. The averaged products of the Green's functions can be now calculated by 
integrating the generating function 



(Ki >m (n,e)) 



■l—m 



11 ml 

We would like to stress that Eqs. 



d* B d*t ( SR (n)s R {n))\s\{n)s A {n)rY{<S> n , $t ). 



(12) 



are valid for an arbitrary hopping matrix 



T. There are only two properties that we require from T at this point: i) T is real and 
symmetric; ii) all the diagonal elements of T are zero. They allow us to rearrange the 
hopping term in the action (flOj) 

^T P9 ($ P! $,) = ^T M Re (s R (p)s R (q) - s A (p)s A (q)) + 



p.q 



p.q 



^2 T pq (x r (p)xr(i) + x* A (p)xA(q)) ■ 



(13) 



v-q 
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The supersymmetric representation of the generating function constructed above is 
exact, but rather complicated, as it involves 8 integration variables (4 real and 4 
Grassmann) associated with each lattice site. Our aim is to demonstrate that in the 
limit e — > one can integrate out 7 out of 8 variables exactly. This can be done by 
changing the variables, which makes it possible to take the limit e — > explicitly in the 
action of the functional integral. A new functional integral representation, obtained in 
such a way, is still exact, but requires an integration over a single real variable on each 
lattice site. 

The first step in realizing this program is to introduce the modulus and the phase 
of the complex variables s R = \sn\e^ H and sa = \sa\g 1 ^ a - Then the non-Grassmann 
part of the hopping term ffl3|) can be written as 

^T pq Re(s* R (p)s R (q) - s* A (p)s A (q)) = ^T pq (cos((j) R (p) - (f) R (q)) \s R (p)\ \s R (q)\ - 

p,q p,q 

cos(<j) A (p) ~ Mq)) Mp)I \s A {q)\) ■ (14) 
Since only the differences between phases appear in the action, we may get rid of the 
phases associated with the nth site by shifting the variables 4> R {p) — > 4> R {p) + ( Pr{ tl ) an d 
4>a(p) —> 0a(p) + 4>A(n) for all p ^ n. This shift does not change the form of the action, 
but it leads to an additional constraint 4> R {n) = 4> A {n) =0. 

The e-dependence of the action ffTO]) suggests that in the limit e — > the main 
contribution to the integrals comes form the regions | s R \ , \s A \ ~ 1 / ^Jt. Therefore it is 
convenient to define two new variables 

s P = e(\s R (p)\ 2 + \s A (p)\ 2 ), v p = \s R (p)\ 2 - \s A (p)\ 2 , (15) 
so that the original variables can be expanded in a power series in e: 

'*y!( i+ S +o <4 ^=M i - e i +o{ ' 2) )' (16) 

where the terms of the order of e 2 and higher play no role in the limit e — > 0. Substituting 
these two expressions into Eq. (fT4|) . one can notice the presence of the large factor 1/e 
in the action, which enable us to integrate over the phases using the stationary phase 



approximation. The details of the calculation are presented in Appendix A and the 
result reads 

J = (2neY N -^Ye lj: ™ TpmVp ^-}— TT - , (17) 

y 1 ^ detS J-l s ' v > 

which is the result of the integration of the non-Grassmann part of the hopping term 
f fl4|) over all phase variables. The matrix B appearing in the above equation is defined 

as 

B pq = -T pq + 5 p q ^ T pr a v a r J—, p, q = 1, . . . , iV; p, q ^ n. (18) 

r V S P 

The discrete variables o~ p — ±1 describe all possible configurations of the stationary 
phase points. Eq.f ll7p represents the leading order in e result for the integral, all higher 
order terms give no contribution to the generating function in the limit e — > 0. 
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The next step is to integrate out the Grassmann variables. However the appearance 
of the quartic terms in the diagonal part of the action ffTO]) prevent us from reaching 
this aim directly. In order to get rid of the quartic terms in the Grassmann variables 



we make a shift of the variables v p : 



v P = v p + X*r(p)Xr(p) + X*a(p)xa(p)- (19) 
This allows us to write the diagonal part of the action ( TlOT) as 



w 



2 



Sdiag = 2^ ( + ^ + S P ) ' (20) 

p V 

where in the last term we retain only the contribution, which gives a non-zero result in 
the limit e — > 0. 

The shift of v p produces new Grassmann terms in Eq. (fl7l) . which together with the 
Grassmann part of the hopping term (fl3il , yields the following integral over Grassmann 
variables: 



where = ^Xr(p) ^Xr{v) ^Xa(p) dXjSp)- The integral now is Gaussian, but its 

action contains the Grassmann variables associated with the nth site, which are not 
integrating out in G. One can eliminate these variables from the integrand by a change 
of the variables similar to the one used in the integration over the phases. Once it is done, 



the integral becomes the standard Gaussian one and the result reads (see Appendix B 
for details) 

G=(-lf- 1 (detB) 2 . (22) 

It is remarkable that the integration over the Grassmann variables leads to the 
appearance of the determinant of exactly the same matrix B as the integration over 
the phases. 

Collecting the results from Eqs.(IT71). ( 1201) and (122]) and taking into account the 
expression for the measure d$d& = —^r t ds dv d<f>Rd<f>AdxRdx*RdxAdxAi we obtain 

Y(s n ,v n ) = J2H r d Vp detB e- E -(^ + ^ +Sp ) +lE - T -^^.(23) 

The variables v p can be easily integrated out, as all the intergals over v p are Gaussian. 
Finally, we introduce the new variables t p = cr PV /s^, such that J2a p Io° ^ = 2 f°° 
then the expression for the generating function integrated over v n reads 



OO tr, ' 



/oo /2,7T r°° dt ~ p \ — 2 — 2 p i 

dv n Y(t n , v n ) = - — TT / p det Be \ ""' > . (24) 

-oo W p^J-oc V^TTWtp 

From Eqs.flO]) and (|12|) it follows that the moments I q (n) are related to Y(t n ) as 

1 f°° 

W> = *mto-*fl L *** nu (25) 
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Eqs. (l24p . (12 5 p along with the definition of B in terms of t p 

B pq = -T pq + dpq^Tpr^, p, q = 1, . . . , JV; p, q ^ n, (26) 

h 

represent the main result of this section. 

The constructed representation for the moments I q (n) is exact and valid for a 
generic hopping matrix T pg . In particular, it can be used to study the original 
Anderson tight-binding model in any dimensions. The main advantage of this approach 
in comparison to the initial representation in terms of the supervectors (110p or the 
supersymmetric non-linear cr-model [18] is that the action depends on a single real- 
valued field t p . Generally, the main obstacle to futher analysis is the presence of the 
functional determinant, which makes the action nonlocal. However, for some models the 
functional determinant can be calculated exactly. One example of this kind is the one- 
dimensional Anderson model. In this case, one can show that the action of the above 
representation can be reduced in the continuous limit to the action of the Liouville 
theory, reproducing the results for I q (n) obtained by the different methods pij[20]. The 
Anderson model on a ci-simplex, considered in the next Section, is another example, for 
which the derived representation turns out to be very useful. 



3. Anderson model on a simplex 

The general approach presented in the previous Section can be applied to the simplex 
model. Since in the derivation of the formula for the generating function (j2l|) we used the 
assumption that T pp = 0, we need to take into account the non-zero diagonal elements 
of the Hamiltonian ([1]) by shifting the energy E — > E — 1/N m Eq. ( |24|) . Once it is done, 
the generating function takes the following form 



,-k £9 -&~ B 



2 



1 \ Jv H L P / 1 .2 

= — n / det Be v • 27 

W p/J-oc V2^wtp 

We are interested in the eigenstates, which are initially degenerate at w = 0. They 
correspond to E = and this is the value of the energy we are going to consider below. 
The functional determinant appearing in the above equation is calculated exactly in 



Appendix C and given by Eq.( IC.15T) . Using this result we obtain 



I q (N) = c N dt n t 2 n g 








-11/ e v 



>2 



N ~ 2 _v f fa** J. 1 fa 



c N = 7= • (28) 

All the sites are equal in the simplex model and therefore we can skip n-dependence in 
the notation for the moments. At the same time we stress the dependence of I q on the 
total number of sites N. 
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Analyzing the action of the functional integral above, we can notice that ^ t q 
plays the role of a collective variable, which determines the behavior of the integral. 
In particular, the condition ^2 q t q = gives a solution of the saddle-point equation in 
the limit w — > 0. This condition can be viewed as a functional integral counterpart 
of Eq.(|3]), which determines the degenerate eigenstates at E — 0. For this reason it 
is convenient to introduce a new variable s by inserting the decomposition of unity 
1 = ds 5(s — (1/wN) ^2„t q ) in the functional integral: 



/oo poo i± poo 

dtnt 2 r 2 \{ / ^ ds8 
-OO £_ J —oo J — oo 




dtp f°° j„ x ( „ \ s iV " 2 e 



1 / x 

C N = i= • 29 

p{0)(q-2)\N{V2^) N w 2 V ' 

Here we extended the lower limit of the integration over t n to — oo by assuming that N 

is even. Now we scale the integration variables t p = sx p and represent the delta function 

by the integral 

wN 




2tt Is 



/oo 
d9e -i0(wN-j: q x q ) (30) 
-oo 



This allows us to perform the integration over all variables x p separately. As a result 
we obtain a very compact representation for I q : 

poo poo 1 

I q (N) = r q / d9e- i6wN / ds \s\ 2q - 3 f N ~\s, 9) g(s, ~ 



J-oo H (2vr)fp(0V(g-2)! 

f(s, 6)= I -^L x- 2 e-^~ s2x2+Wx , g(s, 6) = f° dx x ^-\~^-^^\ (31) 

27T J-oo 



Eq. (l3T!) is exact and it can be used to calculate I q for any finite (even) number of sites 
N. On the other hand, the explicit N dependence makes it possible to analyze the 
integral in the limit iV — > oo and this is the content of the next Section. 

4. Moments of the eigenstates in the thermodynamic limit N — > oo 

In order to find an asymptotic behavior of the double integral (I3T|) as N — > oo we notice 
that the function 

/oo 
ds\s\ 2c >- 3 f N - 1 (s,6)g(s,6) (32) 
-oo 

has a maximum at 9 = 0. Indeed, it follows from the definitions of f(s, 9) and g(s, 9) 
that |/(s,0)| < /(s,0) and \g{s,9)\ < g(s,0), hence 

/oo 
cfeH^I/M)!^ 1 \g(s,9)\<u(0). (33) 
-oo 

It means that in the limit iV — > oo the main contribution to the integral over 9 in Eq. fl3TT) 
originates from 9 — > 0. Therefore it is convenient to introduce the new integration 
variables a = N9 and t = Ns in Eq. (l3T|) : 

W = ^f dae— r dt |t| 2 - 3 f»-> (*°) g (±°). (34) 



N 2 *- 1 J_„ 11 \N'N] \N 7 N 
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Now we need to find the asymptotic expressions for f(et, ea) and g(et, ea) as e — > 0. For 
g it can be done just by changing the integration variable x in Eq.f l3~Tj) by y = x/N: 

t a 
N'N, 

The integral over y can be now explicitly calculated in terms of elementary functions 
t a \ . , r, ,-2q+i ( '•• \ 



9 



N 2q 



-i 



dy y 2q - 2 e~ y 1 +iay + 0(N 



-2^ 



(35) 



9 



N' N 



N 2 q -l j| 



t 



2tJ 



)+0(N- 2 ) 



p=0 



2p)! 



(36) 
(37) 



For /(et, ea) more careful analysis is required, as the function is not analytic at e = 0. 
Its asymptotic expansion at e — > can be written as 

f(et, ea) = f(0, 0) + e lim f' e (et, ea) + 0(e 2 ). (38) 

Calculating the derivative /g(et, ea) using the integral representation for /, which is 
given by Eq. fl3T|) . we obtain 



f(et,ea) = 1 - V2\e\ \t\e~\») 
Applying this expansion with e = l/N to f N_1 (4, -tH we find 



/ 



AT-1 



iV' iV 



-^[i|e-(*) 2 -Vf[«l«rf(|^|) + 0(^- 



(39) 



(40) 



Now we can substitute the results of Eqs.f l36|) and (jiOj) into Eq.f l34|) and take the limit 
N ^ oo 



T q = lim I q (N) = r q 

N— >oa 



* r ^ Fg ( « ) e -^-^|e-(t) 2 -^f H e rf (| f |) 



1 



^-2)!^ \t 



tt / c?2; F„(z) e 



-2iztw—V2\t\e~ 



27r|2i|erf(|z|) 



(41) 



where in the last line we changed the variable a by z = a jit and substituted into the 
expression for r q the density of states calculated in Appendix D 



Since the integral over t looks rather simple, we would like to change the order of 
the integration and to integrate first over t. However, the integral over t diverges at 
t — > 0, so it should be first regularized. One way to do it is to replace 1/ \t\ term by 
1/ integrate over t and then take the limit S — > 0. This procedure yields 



00 F (z) 
-oo (9-2)! 



In 



Az 2 w 2 + 2 (e' z2 + y/rr \z\ exf(\z\) 



? = 2,3,. ...(42) 



Eq.( l4"2"j) and Eq. (157)1 represent the main result of this Section. They allow us to calculate 
the positive integer moments I q at arbitrary value of w. In particular, the fact that the 
moments remain finite in the limit iV — > oo means that the eigenstates are localized at 
any disorder strength w. 

One can check that in limit w — > oo Eq.( l42|) gives I q = 1 for any q, which is in 
agreement with the expectation that at very strong disorder the states are localized at 



Anderson localization on a simplex 



10 




Figure 1. Numerical (symbols) and analytical (solid lines) results for I2 and I4 as a 
function of w. 



single sites of the lattice. On the other hand, the fact that Eq. (I42j) predicts non-trivial 



values for I q at w — > is not so obvious. 



Although in our derivation we assumed that q is integer, one can easily generalize 
the result for I q to non- integer values of q. To this end one needs to calculate the integral 
in Eq.(l35l) assuming that q is non-integer and to replace (q — 2)! in Eq.(l42l by T(q — 1): 



F q = r [ q 



dz 



F q {z) 



T(q 



In 



4z 2 w 2 + 2 



(e- z2 + ^\z\ert(\z\) 



q>l, 



1 1 

q ~ 2'2' 



(43) 



where iFi(a,b,z) is the Kummer confluent hypergeometric function. 

In order to test the validity of Eq.( l4"2"j) and Eq.( l4"3"j) at different values of q and 
w we compare them with the results of numerical simulations. The moments of the 
eigenvectors were calculated by direct diagonalization of 500 x 500 random matrices 
defined by Eq.flTJ). The number of realizations was chosen in such a way, that 5000 
different eigenvectors with eigenvalues close to E = were generated for each value 
of w and q. In Fig. [1] the results for I2 and J4 for various disorder strength w are 
compared with Eq. (j42|) . The numerical and the analytical results for I q as a function 
of q are presented in Fig. [2] for w = 0.01, 3 and 10. Here the analytical predictions were 
obtained from Eq. (|43|) . One can notice that the results of numerical simulations are in 
excellent agreement with the analytical predictions for all values of w and q. 
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q 




Figure 2. Numerical (symbols) and analytical (solid lines) results for I q as a function 
of q for w = 0.01, 3 and 10. 

5. Conclusions 

We have presented a field-theoretical representation for the moments of the eigenstates, 
which involves a single real-valued field variable. The representation is exact and can 
be applied to a generic Anderson model with non-random hopping elements. The 
action of the constructed field theory is non-local, due to the presence of the functional 
determinant. In certain cases, the functional determinant can be calculated exactly. 
Generally, our representation can serve as a starting point for further approximate 
calculations, which might be efficient, for example, in the limit of strong or weak 
disorder. 

We would like to point out that the derived field theory shares some similarities 
with the partition function of the hyperbolic sigma model studied recently in Ref.[21j. 
Indeed, if we restrict all t p in Eq.( l24l) to be positive and replace them by t p = e dp , then 
the action expressed in the new variables 8 P is similar to the action of the the hyperbolic 
sigma model in the horospherical coordinate system. Therefore it might be possible to 
apply the methods from Ref . [21] to prove rigorously the existence of the diffusive phase 
in the three-dimensional Anderson model. 

We have applied the general formula for the moments of the eigenstates to the 
simplex model and derived a very compact representation for the moments. Our analysis 
has shown that the moments remain finite in the thermodynamic limit implying that 
the eigenstates are localized at any strength of disorder. The analytical expression 
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for the moments that we obtained agrees completely with the results of the numerical 
simulations. 

The fact that the eigenstates are always localized might be in conflict with the 
expectation that the presence of the hopping term, which connects any to sites of the 
lattice, must lead to derealization at least at weak disorder. In order to understand, 
why this is not the case, one can consider first the opposite limit of strong disorder. If 
disorder is sufficiently strong, then we expect that the eigenstates are localized. Treating 
the hopping term T as a perturbation, we can notice that it is a rank one matrix, as it 
is given by T = (1/N) \s) (s\, where \s) and (s\ are defined by Eq.f lC.7j) . One can argue 
that the rank one perturbation remains always week independently of the strength of its 
matrix elements. Indeed, one can show that our approach can be used for any rank one 
hopping matrix and the results will be qualitatively the same. We would like to stress 
that the fact, that the hopping elements are non-random is crucial for this conclusion. 

Another way to explain the observed localization of the eigenstates is to notice, that 
in contrast to the standard Anderson model, the eigenstates of the simplex model are 
degenerate at w — 0. The presence of disorder of the strength w leads to the appearance 
of the energy band of the width of the order of w, as it follows from the expression for 
the density of states p ~ 1/w. Hence the bandwidth is always of the same order as the 
disorder strength, implying that disorder is actually strong at any value of w |22j . 



Appendix A. Integration over the phases 

In this Appendix we calculate the integral over the phases 
J = 



2n 



q R 

iJpq 



W d4>n(p) d<p A (p)e i ^^ Tpq ^ cos ^ R ^~^ R ^ 9 p q ~^ 

' 'Vr, 



s p s q ( j £ 
2e V 2 



p _|_ —L 



> 9pq 



y/Spi 



1 



V_p + V_q 



(A.l) 



using the stationary phase approximation, which gives an exact result in the limit e — > 0. 
The stationary phase condition applied to the integral ( lA.ip over 4>r{p) yields 

J2T Pq sm(Mp) ~ Ml))9p q = 0, (A.2) 
q 

1, . . . , N except for p = n. This system of equations 



which should be satisfied for all p 
has obvious solutions 

<Pr{p) - <i>R(q) = n 



if T pq ^ 0. 



(A.3) 



Taking into account that 4>n(n) = 0, we conclude that each variable 4>r{p) has two 
stationary phase values and 7T, which it can take independently from a value of any 
other variable 4>R(q). It is convenient to introduce the new variables 



o-r{p) 



1. 



i*(p) = 

<J>r(p) 







7T 



(A.4) 
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such that the 0#-dependent part of the action can be expanded around the stationary 
phase values as 

9p q cos((j) R (p) - (f) R (q)) = (T R {p)a R {q)y/s^r q x 

/ 1 l I / 11 it \ \ 

(A.5) 



where we neglect all higher order terms in (<f) R (p) — 4> R {q)), as they give no contribution 
in the limit e — > 0. 

Repeating the same steps for the integral over 4>a(p), we obtain 

9pqCOs((j)A(p) ~ <j>A(q)) = 0-A(p)v A (q)y/Sp~Sq~ X 

s-5<^-^>> 2 -ife+3)> < A '«) 

Comparing these two results, we notice that the singular 1/e contribution in the action is 
cancelled out iff c R {p)a R {q) = cr y i(p)(j^(g). Taking into account that o~ R {n) = cta(^) = 1 
by definition, we conclude that a R {p) = a A (p) = a p for all p. 

For a fixed configuration of {c p }, the integral over <fi R {p) is given now by the 
Gaussian integral in the leading order in e: 

I" T]dMp)e-^-^^ iMp) - Mq))2 = ( 6 ~ tJ ^} . (A.7) 

The matrix Q is a (N — 1) x (N — 1) matrix defined as 

Qpg = o- p o- qy /Sp~s^ ^-T pq + 8 pq ^2 T prO-pO-r ^j-^j , p, q = 1, • • • , N; p, q ^ n. (A.8) 

One can notice that Q can be factorized as 

Q = (A.9) 
where S p(? = o- py /s~p~5 pq and 

B pq = -T pq + 5 pq ^2T pr a p cr r , — , p,q = 1, . . . , N; p,q n. (A. 10) 



p 



The factorization enables us to express the determinant of Q as 

det Q = det B J[ s p . (A. 11) 

The integral over <Pa{p) can be evaluated in a similar way: 

/ -7T / \N~ 1 

r rr^ A ( P )e^^^^v^(^(p)-^w) 2 = ( e * T v^f[ — (A12) 

Collecting the results from Eqs. (lA.5j ). (1A.7j) and ( 1A.12|) we find 



J = (2*e)W-V PqWP ^^B II ~ p {A.13) 

{o-p} p^n 
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We assume that the system of equations (1A.2j) has only trivial solutions given by 
Eq. (IA.3j) . which is the case for a generic matrix T pq . If for a specific choice of T pq a 
non-trivial solution of ( 1A.2I) does exist, its contribution to the integral should be also 
taken into account in a similar way. 



Appendix B. Integration over the Grassmann variables 

The integral over the Grassmann variables ( 12T1) can be factorized as follows: 

p^n 



Cpq Tpq Spq ^ ^ Tp r O~pO' r 



(B.r 



This integral is Gaussian, but the quadratic form Yli Pq CpqX*pXp contains variables Xn 
and x*ni which are not integrated out. In order to eliminate these two variables from 
the integrand, one can perform the change of the variables 

Is 

Xp = r] p + DpXn, X* P = V*p + D pXl, Dp = o p a n J^. 

V &n 

In terms of the new variables the quadratic form reads 

Y C PlX*pXq = Y GpqVpVq + Y C pq D qVp + Y C P^lp Xn + 



(B.2) 



p,q 



I"! 



p,q 



X n I ^ ] GpgDpTjq + ^ ] CpnVp I + I ^ CpqD p D q + 2 ^ ] C pn D p + C nn 

VnVn, (B.3) 

V p,q p / \ p,q p / 

where the symbol ^2 m eans that the terms with p = n and q = n must be excluded 
from the sum. 

Now we show that all the terms, except for the first one, on the right hand side of 
the equation above are equal to zero. Indeed, it follows from the definition of C pq , that 
for p ^ n 



^ GpqDq ~\~ Cpn ^ ^ | Tpq &pq ^ Tp r O~pO~ r 



O q O n 




+ T 



pn 



■L pqO~q®~n ~ / -L pr^r^n \j ~ "T -t pn 



(B.4) 



Hence the second and the third terms in Eq. (lB.3j) vanish. Similarly for the last term in 
Eq.( 1B.3l) we obtain 

( 

^ CpqDpDq ^ Cpq | T pq 5pq ^ Tp r O pO ~ r ^ 



p,q 





pqOpO q - 



' SpSq 



prO~pO~ r 



SpSq 



— T m a p (j n J — . (B.5) 



p,q 



p r 
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On the other hand the sum of the remaining two terms is equal to 
/ / , — — 

2 ^ ^ C.p n Dp -\- C nn 2 ^ ^ Tp n OpO~ n \l — -\- T nn ^ ^ T rn o~ r o~ n ^ 




y^T pn o-pQ- n J — , (B.6) 

V 

so that we conclude that the last term in Eq.( 1B.3l) equals to zero. Thus we find from 
Eq.flHU) 



Yl C piX*p*p = C Pi r ll r li = ~ B v<fh^ (B.7) 
p,<? p,q p,q 

where B pq is defined by Eq. (fl8|) . The integral in (IB.lj) is now the standard Gaussian 

one and we immediately obtain the result 

G={-l) N -\detB)\ (B.8) 
Appendix C. Calculation of the functional determinant 

The matrix B pq defined by Eq. (T26|) is the (N — 1) x (N — 1) matrix, which can be 
obtained from the N x N matrix B 

B pq = -T pq + 5 pq Y TprJ 1 , P,q=l,...,N, (C.l) 

by eliminating all the elements in nth row and in the nth column. The matrix B is 
singular, as the vector 

N \ -V2 



Z 



Yi) (tut 2 ,...,t N f (C2) 

vP=l 



is the normalised eigenvector of B corresponding to the zero eigenvalue. 

Using the first order perturbation theory for the eigenvalues, it is easy to show that 

det B n = {-l) p+c >z;z q detS, (C.3) 

where B pq is an (N — 1) x (N — 1) matrix, obtained from B by eliminating all the 
elements in pth row and in the gth column and deti? is the product of all non-zero 
eigenvalues of B. This formula is valid for an arbitrary singular matrix B, which has 
the non-degenerate zero eigenvalue with an eigenvector z. In particular, it implies that 

det 5 = \z n \ 2 detB (C.4) 

On the other hand, the product of all non-zero eigenvalues of B can be calculated as 

detB = lim-detCB + eI), (C.5) 

where I is the N x N identity matrix. This result simply follows from the fact that 
adding of el to B shifts all the eigenvalues of B, including the zero one, by e. 
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For the simplex model the matrix B reads 



B 



Nt v <Er=l l r N ' P ~ " 



pq 



77 > P^Q 



so that one can represent it using the Dirac notation as follows 
1 



B 



N 



\s) (s\ + aD), 
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(C.6) 



(C.7) 



A' 



(s\ = (1, !,...,!), |s) = (l,l,...,l) T , a = J>, D = die,g(t^,t^,...,t N 1 ). (Ci 



r=l 



According to Eq. (lC.5p we need to compute 



(C.9) 
(CIO) 

where e stands for Ne. The first determinant on the right hand side of Eq. (1C.9j) can be 
found easily, as the matrix el + aD is diagonal: 



det (B + el) = — w det(el + aD) det (J - G \s) (s|), 
G = {el + aD)-\ 



N , v 

det(eI + aD) = JJ (- + ej . 

p=i ^ p ' 



(en; 



In order to calculate the second determinant on the right hand side of Eq.( lC.9j) we notice 
that 

lndet(J - G Is) (s|) = tr ln(/ - G Is) (s|) = - V - tr (G Is) (s\) n = 

^— ' n 

n>l 

- E l^ s \ G ° i s » n = ln ( x - ^ G ° I s )) = ln ( 1 - E (f + ") 

n>l \ p=l ^ p ' 



(C.12) 



From Eqs.flDU, (ICTTj) and (jUl2]l we find that 



det (5 + eJ 



P =i 



a 



Then it follows from Eq. (lC.5p that 



— — ' , 1 1 
det-£> = hm - 



N / N -l\ AT 

>-e ?« n 



p=l 

N 



N-2 



e^o e N N 



p=l 



=1 v*p 



( 7, 



+ e 



(C.13) 



p=i 



(C.14) 



r=l \r=l 

Substituting this result into Eq. (lC.4|) and taking into account Eq.( lC.2j) we obtain finally 
the expression for the determinant of the matrix B: 



N-2 



det B 



, 2 n \ »-* N 



(C.15) 



p=l 
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Appendix D. Calculation of the density of states 

The averaged density of states is defined as 

1 N 1 
p=i 

where iO is an infinitesimal imaginary shift in energy and the Green's function G is given 
by 

For the simplex model the Green's function reads 

G(S)=tT(G^-\s)(s\)-\ (D.3) 

where (Gq = 5 pq (£ — H PP )N and \s) and (s\ are defined by Eq. (IC.7j) . Since Go is 
diagonal and \s) (s\ is an operator of rank one, the inverse operator in Eq.( lD.3l) can be 
calculated explicitly: 

G{£) = tr(/ - Go \s) (s\)G = tr ^(G \s) (s\) k G = 

k>0 

tr (Go + T((s\ Go |s)) fc - 1 G \s) (s\ G J = trG + - ( ''^ . (D.4) 

The first term here corresponds to the diagonal part of the Hamiltonian and its 
contribution to the density of states is equal to P(E), where P(H PP ) is the probability 
distribution of the diagonal elements. One can show that the contribution of the second 
term vanishes in the limit iV — > oo due to the presence of Gq in the numerator. As a 
result we obtain 

p(0) = P(0) = -=^. (D.5) 

V Z1TW 
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